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We present a first-order symmetric hyperbolic system in the Ashtekar formulation of general relativ- 
ity for vacuum spacetime. We add terms from the constraint equations to the evolution equations 
with appropriate combinations, which is the same technique used by Iriondo, Leguizamon and Reula 
( Phys. Rev. Lett. 79, 4732 (1997) ). However our system is different from theirs in the points 
that we primarily use Hermiticity of a characteristic matrix of the system to characterize our system 
symmetric, discuss the consistency of this system with reality condition, and show the characteristic 
speeds of the system. 

PACS numbers: 04.20.Cv, 04. 20. Ex, 04.20.Fy 



I. INTRODUCTION 

Hyperbolic formulation of the Einstein equation is the 
one of the main research areas in general relativity |Q. 
This formulation is used in the proof of the existence, 
uniqueness and stability (well-posedness) of the solutions 
of the Einstein equation by analytical methods H . So far, 
several first order hyperbolic formulations are proposed; 
some of them are flux conservativ e f| , symmetrizable [Q , 
or symmetric hyperbolic system []5h7|. The recent inter- 
est in hyperbolic formulation arises from its application 
to numerical relativity. One of the expected advantages 
is the existence of the characteristic speeds of the system, 
with which we may treat the numerical boundary with 
appropriate condition. Some numerical tests have been 
reported along this direction |^-|Tc[|. 

Recently, Iriondo, Leguizamon and Reula (ILR) Q 
discuss a symmetric hyperbolic system in the Ashtekar 
formulation |l2| of general relativity. Ashtekar 's formu- 
lation has many advantages in the treatment of gravity. 
By using his special pair of variables, the constraint equa- 
tions which appear in the theory become low-order poly- 
nomials, and the theory has the correct form for gauge 
theoretical features. These suggest possibilities for treat- 
ing a quantum description of gravity nonperturbatively. 
The classical applications of the Ashtekar's formulation 
has also been discussed by several authors. For example, 
we [Q discussed the reality conditions for metric and 
triad and proposed new set of variables from the point 
of Lorentzian dynamics. We fig ] also showed an example 
of passing degenerate point in 3-space by loosing reality 
condition locally. 

In this Letter, we present a new symmetric hyperbolic 
system in the Ashtekar's formulation for Lorentzian vac- 
uum spacetime. ILR 11 says that they construct a 



discussion is not clear in the the following three points. 
First, they used anti- Hermiticity of the pricipal symbol 
for defining their system symmetric. We, however, think 
that this does not derive He rmit icity of the characteris- 
tic matrix [A below in eq. (3.1), B in pd{|], since they 
do not define their vector k a explicitly. We rather use 
Hermiticity of the characteristic matrix primarily to con- 
struct a symmetric hyperbolic system. Second, they did 
not mention the consistency of their formulation with 
the reality conditions which are crucial in the study of 
the Lorentzian dynamics in the Ashtekar variables [fl9|| . 
Third, they did not discuss the characteristic structure 
of the system, which should be shown in the normal hy- 
perbolic formulations. Our discussion covers these two 
matters. 

The construction of this paper is as follows. After giv- 
ing a brief review of Ashtekar's variables and reality con- 
ditions in §2, we present our formulation in §3. The dis- 
cussion of characteristic speed and summary are in §4. 



II. ASHTEKAR'S FORMULATION 

The key feature of Ashtekar's formulation of general 
relativity |l2| is the introduction of a self-dual connection 
as one of the basic dynamical variables. Let us writef] 
the metric g^ using the tetrad, e^, and define its in- 



E i, by 



^ivu and 



We 



symmetric hyperbolic system. However, we think their 



* We use /i, v — 0, • ■ ■ , 3 and i,j — l,---,3 as spacetime 
indices, while L,J = (0), • ■ • , (3) and a,b — (1), • ■ ■ , (3) are 
50(1,3), SO(3) indices respectively. We raise and lower 
■ ■ • by g^ v and g^ v (Lorentzian metric); 7, J, • • ■ by rj IJ = 
diag( — 1, 1, 1, 1) and rju; ■ ■ ■ by 7 lJ and jij (3-metric) . We 



use volume forms e a bc, e a bc£ a 



3!. 
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define S0(3,C) self-dual and anti self-dual connections 
± A!f l := w" a ^(z/2)e a bc w^ c , where uj 1 ^ is a spin connection 
1-form (Ricci connection), uoj/ := E Iv V M e^. Ashtekar's 
plan is to use only a self-dual part of the connection + A C ^ 
and to use its spatial part + Af as a dynamical variable. 
Hereafter, we simply denote + A a ^ as A^. 

The lapse function, N, and shift vector, N l , are ex- 
pressed as Eft = (1/N, -N*/N). This allows us to think 
of Eq as a normal vector field to £ spanned by the con- 
dition t = x° =const., which plays the same role as that 
of ADM. Ashtekar treated the set (Af, El) as basic dy- 
namical variables, where E l a is an inverse of the densitized 
triad defined by E\ := eE l a where e := det is a density. 
This pair forms the canonical set. 

In the case of pure gravitational spacetime, the Hilbert 
action takes the form 

S = [ dM(d t At)K + (t/2)NElElF^ ab c ~ANdetE 



-N'F^Ei + A^ViEi] 



(2.1) 



where N := e X N, A is the cosmological constant, 

V t E l a := d, t El - ie ab c A b E l c , and det£ is defined to be 
detE = (l/6)e o6c e iife ^^'4 fc , where e ijk := e abc e?e<>e% 

and e.ij k ■= e~ 1 e ijk |. 

Varying the action with respect to the non-dynamical 
variables N, N l and Aq yields the constraint equations, 



C H := (i/2) e"" c E l a El Ff- - Adet£ a 0, 
Cm* : / « 0, 
C Ga := 2?i£* w 0, 



(2.2) 
(2.3) 
(2.4) 



where := (dA a )^ - (i/2)e Q &c (^ b A A c )^ is the cur- 
vature 2-form. 

The equations of motion for the dynamical variables 
(A? and E l a ) are 



Notice that in general the metric ( |2.7| ) is not real. To 
ensure the metric is real-valued, we need to impose real 
lapse and shift vectors together with two conditions (met- 
ric reality condition); 



Jm(E i a E ja ) = 0, 
Re(e abc E*E { b l V k Eiy) = 0, 



(2.8) 
(2.9) 



where the latter comes from the secondary condition of 
reality of the metric lm{d t (Ei& a )} = p), and we 
assume detE > (see pi|). 

For later convenience, we also prepare stronger reality 
conditions. These conditions are 



Im(El) = 
and Im(d t E l a ) = 0, 



(2.10) 
(2.11) 



and we call them the "primary triad reality condition" 
and the "secondary triad reality condition" , respectively. 
Using the equations of motion of JEJ , the gauge constraint 



4|), the metric r eality conditions (2 



.91) and the 

primary condition (2.10), we see that (2.11) is equivalent 
to|| 

Re(A%) = di(N)E la + {l/2e)e b NE 3a d j E l b + N l Re(A?), 



or with un-densitized variables, 



Re(A%) = di(N)E ia + N^eiAf). 



(2.12) 



(2.13) 



From this expression we see that the second triad re- 
ality condition restricts the th ree c omponents of "triad 
lapse" vector Aq. Therefore (2.12) is not a restriction 
on the dynamical variables (Af and E l a ) but on the slic- 
ing, which we should impose on each hypersurface. Thus 
the second triad reality condition does not restrict the 
dynamical variables any further than the second metric 
condition does. 



d t A« 
dtEi 



b NEiF c 

c rZ b %3 
cb 



(2.5) 



III. HYPERBOLIC FORMULATION 



W^NEiEl) + W^ES) + iAleJEl (2.6) 



where tyX? := djXl 



X? = 0. 



Kb-^XiyoTXv , 
In order to construct metric variables from the vari- 
ables (Af,El,N,N l ), we first prepare tetrad Ej as 



E 



,, - {l/eN,-N*/eN) and£# = (0,E*Je). Using them, 



we obtain metric such that 



f v := ElE v jT] IJ . 



(2.7) 



6, €xyz — 1) f: 



We start from defining hyperbolic system following 
Friedrichs (T^| , which is first applied in general relativity 
by Fischer and Marsden ||. That is, we say that the 
system is first-order (quasi-linear) hyperbolic if a certain 
pair of variables Ui form a linear system as 



dtiit = A^^diUj + Bi(u), 



(3.1) 



where A is a characteristic matrix-valued function, of 
which eigenvalues are all real, and B is a function. We 
further define that the system is symmetric when A is a 
Hermitian matrix 16^1^]. 

The symmetric system gives us the energy integral in- 
equalities, which are the primary tools for analizing well- 
posedness of the system. As was discussed by Geroch 
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p8| , most physical systems are expressed as symmetric 
hyperbolic systems. 

Ashtekar's formulation itself is in the first-order form 



in the sence of (3.1), but not a symmetric hyperbolic 
form. 

We start from writing the principal part of the 
Ashtekar's evolution equations as 







Af_ 





A I bi 
A a j 

° ij u bi 



B ? 

ab 
-da j 



(3.2) 



where = means that we extracted only the terms which 

appear in the principal part of the system. The system 
is symmetric hyperbolic if 

= A labij - A lbaji , (3.3) 

= D labij -D lbaij , (3.4) 

Q = B labij _ Qlbaji ( 3 5 ) 

where bar denotes complex conjugate."'" 

We first prepare the constraints (pT^)-(2.4) as 

C H = ie ab c ElEid^ = le^E^A)) 

= -le^El&MA)), (3.6) 

C Mk = -F^El S ~{d k A a 3 - d,A a k )Ei 

= [S l k El+5iE l b ](d l A b ), (3.7) 

C Ga = ViE l a S 8iE l a . (3.8) 



We apply the same technique with ILR to modify the 
equation of motion of E a and A\ by adding the con- 
straints which weakly produce Ch = 0, Cuk = 0, and 
Caa = 0. With a parametrization for triad lapse Aq with 
T and S as 



a i j%&T l i a i b dt& b + s l i aj b d l j& 



(3.9) 



we write the principal parts of (2.5) and ( |2.6[ ) as 

3 t E a = -W^NEiEl) + 2V 3 {N^E*}) 
+iA b e ab c E l c + P\ b C b G 
S -ie^NidMWb ~ i^aNE^El) 
+V j {WE i a ) - V 3 {N*El) + P\ b d 3 El 



rb 



NE[5) 



[-ie\N8\Ei- i, 
+N L S)S b a ~N%S b a +P\ b S l A(d l El), 



(3.10) 



d t A? = -ie ab NEiFZ + WF? i 



J We think that the reader will not confuse A labij and B labi] 



with matrix A and B in (3.1) 



+V. t A a + eANel + Q^C H + R/ a C M] 

= -ie ab c NE 3 b (d l A'j - djAt) + N j (djAi - d.A^) 

+T\ a 3 - b d l & b + s';'\<hA ! ] - Q a k ie b cd E j c E l d {diA b ) 
+R i ka [-S{Ei+SiE{]d l A b 
= [+ie a b c NEl5\ - ie a b c NE l c 8i + N l S a b S{ - N^S\ 

,ql aj 

+D i b~ 

l|al,n " 5i , (3.11) 



where P, Q and R are paramet ers an d will be fixed later. 
Note that we truncated Aq in ( 3.1C| ), while it remains in 
( |3.11 ), since only the derivative of Aq effects the principal 
part of the system. From these two equations, we get 



j^abij 
j^labij 



-ie bca N^E l c - ie cba NE l cl lj 

+N l Y j S ab - N i Y j S ab + P iab j l i } 
0, 

^jlabij rj-iliajb 

D labij = +ie »bc^p cl li _ i^bcy^ji 

+N l 5 ab j ji - N J S ab j u + S lia]b 
-iQ al € bcd E j r E l H - R Ua E 3b + R ija E lb . 



(3.12) 
(3.13) 
(3.14) 



(3.15) 



The condition (3.3) is written as 



= -ie^NY^El ~ ie acb NY l E c 
-2ie cba Nlm{E l c )Y j - N i -y lj 5 ab 
+N j -f lt 5 ba + . 



piab^l] _ p>„ 



ba_li 



(3.16) 



Because the third term in the right-hand-side cannot be 
clliminatcd using P, we assu me th e tria d rea lity condition 
Im(E l c ) = hereafter. Then Q) and (§1J) become 



(3.17) 



= J^ljbia 

= ie abc NEij H + ie bac NE^ - N j 5 ab j H + A^ b< V 3 

+S liajb _ gljbia _ iQriJcdp^ty _ iQ b i e acd E' c E d 

-R Ua E jb + R l ^ a E lb + W lb E ia - R jtb E la . (3.18) 



The third and forth term in ( |3.18 ) cannot be elliminated 
using Q or R, so S Uajb = 7 ' l <F r A" is determined. Thus 
S and T l i Ua = [eq. ( |3.17D ] decides the form of the triad 
lapse as 



Aq = A°jN 3 + non-dynamical terms 



(3.19) 



in result. In or der to be consistent with the triad lapse 
condition ( 2.13 ), we need to specify the lapse as diN = 
0. This lapse condition is also supported by the fact 
that if we do not assume diN = 0, then the secondary 
triad reality condition ( 2.12| ) makes the system second 
order. ILR does not discuss consistency of the system 



3 



with reality condition (especially with secondary reality 
condition). However, since ILR assume Aq = A^N^ , we 
think that ILR also needs to impose similar restricted 
lapse condition in order to preserve reality of the system. 

The rest of our effort is finished when we speci fy pa - 
rameters P, Q and R. P is given by decomposing ( 3.16| ) 
into real/complex parts; 







-N i -y lj 8 ab + N j j H 6 ba 



+Re(P) w V J 



Re(P) ifca 7 i ' 



- e bca NY j Ei - e acb Nj H Ei 



+Im(P) Mb 7 y + Im(P) ifc V l 



(3.20) 



(3.21) 



By multiplying 7^ in these two and taking symmetric 
and anti-symmetric operation to the index ab, we obtain 

piab = N t S ab + j7y e a6c_gi < ^2) 

For Q and R, we found that a combination of the choice 
qai = e -2jy^i Q (3 23) 

R ila = i e -2jy- e acd^i^ (3 24 ) 

satisfies the condition ( 3.18|) . 



IV. DISCUSSION 

In summary, by adding constraint terms with appro- 
priate coefficients, we succeed to construct a symmetric 
hyperbolic formulation for the Ashtekar's system. This 
formulation is consistent with secondary triad reality con- 
dition, which requires to impose a constant lapse function 
for the evolving system. 

The characteristic speeds of this system are give n b y 
finding eigenvalues of the characteristic matrix A of (3.1). 
Since A is a Hermitian, eigenvalues of A are all real. Then 
it is again clear that this system is symmetric hyperbolic 
Actually the eigenvalues of the 18 x 



18 matrix A 1 for 



ar -direction are: N l (multiplicity = 6), N l ± y/j u N (5 
each), and N l ± 3y // j^N (1 each), where we do not take 
the sum in j u here. These speeds are independent from 
the way of taking a triad. We omit to show the related 
eigen- vectors because of saving space. 

As we denoted in §3, our formulation requires triad 
reality condition. In order to make the system first or- 
der, the lapse function is assumed to be constant. Shift 



vectors and triad lapse Aq should have a relation ( p,19[ ) . 
This can be interpreted that shift is free and triad lapse 
is determined. This gauge restriction soun ds ti ght, but 
this arises from our general assumption of ( |3 . 9[ ) . There 
might be a possibility to improve the situation by renor- 
malizing shift and triad lapse terms into left-hand-side of 
equations of motion like the case of GR Q] . Or this might 
be because our system is constituted by Ashtekar's orig- 
inal variables. We are now trying to release this gauge 



restriction and/or to simplify the characteristic speeds by 
other gauge possiblities and also by introducing new dy- 
namical variables. This effort will be reported elsewhere. 

This work (HS) was partially supported by NSF PHYS 
96-00049, 96-00507, and NASA NCCS 5-153. 
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